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1 Introduction 

The separation of variables for solving by quadratures the Hamilton-Jacobi (HJ) equations of 
related Liouville integrable dynamic systems with quadratic in momenta first integrals has a 
long history as a part of analytical mechanics. There are some mile stones of that theory. First, 
in 1891 Stackel initiated a program of classification of separable systems presenting conditions 
for separability of the HJ equations in orthogonal coordinates [l]-[3]. Then, in 1904 Levi-Civita 
found a test for the separability of a Hamiltonian dynamics in a given system of canonical 
coordinates [4]. The next was Eisenhart [5]- [7], who in 1934 inserted a separability theory in 
the context of Riemannian geometry, making it coordinate free and introducing the crucial 
objects of the theory, i.e. Killing tensors. This approach was then developed by Woodhouse 
[8], Klanins [9]-[10] and others. Finally, in 1992, Benenti [11] -[13] constructed a particular but 
very important subclass of separable systems, based on the so called special conformal Killing 
tensors. 

The first constructive theory of separated coordinates for dynamical systems was made by 
Sklyanin [14]. He adopted the method of soliton theory, i.e. the Lax representation, to systematic 
derivation of separated coordinates. In that approach involutive functions appear as coefficients 
of characteristic equation (spectral curve) of Lax matrix. The method was successfully aplied 
to separation of variables for many integrable systems [14]- [21]. 
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Recently, a modern geometric theory of separability on bi-Poisson manifolds was constructed 
[22]- [29], related to the so-called Gel'fand-Zakharevich (GZ) bi-Hamiltonian systems [30], [31]. 
Obviously, it contains as a special case Liouville integrable systems with all constants of motion 
being quadratic in momenta functions. Indeed, Ibort et.al. [32] proved that the Benenti class of 
systems can be lifted to the GZ bi-Hamiltonian form. 

In the following paper we construct in a systematic way a separability theory of all Liouville 
integrable systems on Riemannian manifolds, which are of the GZ type, including as a special 
case the Benenti class of systems. Actually, first we construct a quasi-bi-Hamiltonian theory 
on the so called ojN manifolds and then we lift it to the related GZ bi-Hamiltonian systems 
on bi-Poisson manifolds. What is important, infinitely many classes of separable systems are 
constructed from appropriate deformations of the Benenti class of systems. In that sense we 
demonstrate the crucial role of this particular class of systems in the separability theory of 
dynamic systems on Riemannian manifolds. 

The organization of the paper is as follows. In Section 2 we sketch the geometric separability 
theory of GZ bi-Hamiltonian systems, which was recently constructed and is a main tool used 
in the paper. In Section 3 we recall the basic facts about separable dynamics on Riemannian 
manifolds. Section 4 deals with a special case of separable systems, i.e. the so-called Benenti 
systems. We review this class of systems systematically as it plays a crucial role in a separability 
theory of Gel'fand-Zakharevich systems on Riemannian manifolds and is of a special importance 
for the theory developed in this paper. In Section 5, we construct the simplest new classes of 
separable systems being the so-called 1-hole deformations of the Benenti class. For this example 
we explain the main ideas of our approach as well as the methods of systematic construction of 
separable potentials and (quasi)-bi-Hamiltonian representations. In Section 6, we develop the 
approach to the case of arbitrary k-hole deformations of the Benenti class of systems, construct- 
ing a complete theory of separable GZ systems on Riemannian manifolds. Finally, in Section 7, 
we illustrate our theory by a few simple examples. 

2 Gel'fand-Zakharevich bi-Hamiltonian systems and their sep- 
arability on ujN manifolds 

2.1 Basic definitions 

In this Section we sketch the basic concepts of the modern geometric separability theory to 
make other parts of the paper more clear to the reader. Let us first remind few basic facts 
from the Poisson geometry. Given a manifold Ai, a Poisson operator tt on M. is a mapping 
7r : T*M — > TM that is fibre-preserving (i.e. tt\t*m ■ T*M — > T X M for any x G M) and such 
that the induced bracket on the space C°°(A / J) of all smooth real- valued functions on M 

{., . } n : C°°{M) x C°°(M) -> C°°(M) , {F, G} n = (dF, tt dG) (2.1) 

(where (.,.} is the dual map between TM and T*M). It is skew-symmetric, satisfies Jacobi 
identityand the Leibniz rule and the symbol d denotes the operator of exterior derivative. The 
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second order contravariant tensor field tt can always be interpreted as a bivector, ir € A 2 (Ai) 
and in a given coordinate system (x 1 , . . . , x rn ) on M. we have 

-E^A^. (2.2) 

The function C : M — > R is called the Casimir function of the Poisson operator 7r if for an 
arbitrary function F : M — ► R we have {F, C} n = (or, equivalently, if ndC = 0). A linear 
combination ir^ = m — ^7ro (£ £ R) of two Poisson operators 7ro and m is called a Poisson pencil 
if the operator ir^ is Poissonian for any value of the parameter £. In this case we say that ttq 
and 7Ti are compatible. Given a Poisson pencil tt^ = tt\ — $,ttq we can often construct a sequence 
of vector fields X; L on M that have a twofold Hamiltonian form (the so-called bi-Hamiltonian 
chain) 

Xi = Tt x dhi = ir dh i+1 (2.3) 

where hi : M —>■ R are called the Hamiltonians of the chain (2.3) and where i is some discrete 
index. This sequence of vector fields may or may not truncate (depending on the existence of 
Casimir functions). 

Let us consider a bi-Poisson manifold (M, ttq, tt±) of dimM = 2n + m where ttq, tt\ is a pair 
of compatible Poisson tensors of rank 2n. Moreover we assume that the Poisson pencil admits 
m, polynomial with respect to the pencil parameter £, Casimir functions of the form 

rij 

h (j)^ )= J2h?e j -\ 3 = l,..,m, (2.4) 

i=0 

(7) 

such that m + ... + n m = n and h\ are functionally independent. The collection of n bi- 
Hamiltonian vector fields 

7r ? d/j(%) = X\ j) = mdh^ = irodh^!, i = 1, -, rij, j = 1, -, m, (2.5) 

is called the Gel'fand-Zakharevich system of the bi-Poisson manifold M. Notice that each chain 

(i) 

starts from a Casimir of ttq and terminates with a Casimir of tt\. Moreover all h) pairwise 
commute with respect to both Poisson structures 

= (dh^^dh^) = (dh^^dh?) = 0. (2.6) 
2.2 Geometric separability theory 

(i) 

The necessary condition for separability of Hamiltonian functions hf is a projectibility of the 
Poisson pencil ir^ onto a symplectic leaf A/" of 7ro, that is a 2n-dimesional submanifold defined by 
fixed values of Casimir functions of tto : h^ = c±, h^ = c m . Thus, M is a submanifold of a 
codimension mm. M. The projection is done through the particular realization of the Marsden- 
Ratiu scheme [33], i.e. along the appropriate transversal distribution. Let Zj, % = l,...,m be 
some vector fields transversal to M, spanning an involutive (integrable) distribution Z in M of 
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a constant dimension m (that is a smooth collection of m-dimensional subspaces Z x C T X A4 at 
every point x in A4). The word 'transversal' means here that no vector field Zi is at any point 
tangent to the submanifold M passing through this point. Hence, the tangent bundle TM. splits 
into a direct sum 

TM=TN®Z (2.7) 

(which means that at any point x in Ai we have T x A4 = T X N © Z x with s such that x G J\f ) 
and so does its dual 

T*M = T*N®Z\ (2.8) 

where T*M is the annihilator of Z and Z* is the annihilator of TJ\f . It means that if a is a one 
form in T*N then a(Zi) = for all i = 1, . . . , m and if (3 is a one- form in then /3 vanishes 
on all vector fields tangent to M. Moreover, we assume that the vector fields Zi which span Z 
are chosen in such a way that dh$\ j = 1, m is a basis in Z* that is dual to the basis Zi of 
the distribution Z, 

(dh\i\z^ = Z l {h { i ) ) = 5i v (2.9) 

(it is no restriction since for any distribution Z transversal to M we can choose its basis so 
that (2.9) is satisfied). The 'orthogonality' condition (2.9) together with the involutivity of the 
distribution indicates the relation 

[Zi,Zj]=0, (2.10) 
where [., .] stands for the commutator of vector fields. 

Theorem 1 Poisson tensor 7Ti is projectible onto a symplectic leaf N o/ttq along a transversal 
integrable distribution Z if the normalized vector fields Z; L locally generating Z are symmetries 
of no (Lz^o = 0) and satisfy 

L Zt 7r 1 = Y,Y t U) AZ„ (2.11a) 

j 

where 

= tto d(Zi(hP)) = [Z^ x dhf\ = [Z^X^] (2.11b) 

and Lz means a Lie derivative in the direction of Z . The projection of tt\ onto M along Z is 
equivalent to the projection of the deformed Poisson tensor 

tt 1d = TTi - x ? A z ii x ? = 7ri^o j) = Kodh^ (2.12) 
j 

onto its symplectic leaf M. Moreover, tti and ttid are compatible. 

The proof as well as the details of the whole construction the reader can find in [29], [34] 
and [35]. Notice that 

L Zi lTlD = 0, 

so Zi are also symmetries of it\D- 
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Let us denote the projections of ttq, tt\ by 9$, 6± and restrictions of (h\ , hh^ )\m to M by 
(H^~\ H^) := (Hi, H n ). Notice that the projection of the Poisson pencil 7T£ is a Poisson 
pencil 0£ on and in a generic case both Poisson tensors 9q,9\ are nondegenerate. Hence, TV 
is a bi-symplectic manifold as is endowed with two symplectic forms uq,ui defined by 

{F,G} ei =Ui(X F ,X G ), X F = 9 dF, i = 0,l. (2.13) 

This simply means that ojq = 9q 1 and coi = uj 8iujq. Indeed 

u (X F ,X G ) = <u X F ,X G >=<uj 9 dF,9 dG > (2.14a) 
= <dF,9 dG>={F,G}g , 

uji(X f ,X g ) = <uq0iu o X f ,X g >=<uoX f ,6iujqX g > (2.14b) 
= <dF,e 1 dG>={F,G} 6l . 

Moreover, one can construct the tensor field N := ^i^ 1 = Oiluo, of type (1, 1), called a recursion 
operator of TV and its dual N* = ujqOi .Notice that 

N9 = 9i, N*uj = uj 1 . (2.15) 

The important property of N is that its Nijenhuis torsion 

T(N)(X, Y) := [NX, NY] - N([NX, Y] + [X, NY] - N[X, Y]) (2.16) 

vanishes as a consequence of the compatibility between #o an d 9\ and hence implies that uji 
is closed [36]. Such manifolds are known as the so-called ujN manifolds. The generic case 
means that 2n-dimensional toN manifold is endowed with a recursion operator N which has 
at every point n distinct double eigenvalues Ai,...,A n , which are functionally independent on 
J\f. Choosing Aj as the canonical position coordinates, we can always supplement a set of local 
coordinates (Aj ;/ Uj) on TV by the canonically conjugate momenta ^j. 

Definition 2 A set of local coordinates (A^/uJ on loN manifold M is called a set of Darboux- 
Nijenhuis (DN) coordinates if they are canonical with respect to 9q and diagonalize the recursion 
operator with diagonal elements beeing its eigenvalues. 

It means that in the (A, /Lt) coordinates 

*=U o)- o> -(»:)• <2i7) 

where A n = diag(Ai, A n ), and their differentials span the T*M which is an eigenspace of N* 
(the adjoint of N), as 

N*d\ = XidXi, N*dfi { = Xidfii, i = 1, n. (2.18) 
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As well known, n functionally independent Hamiltonian functions Hi,i = 1, ...,n are said to 
be separable in the canonical coordinates (A, fi) if there are n relations, called the separation 
conditions (Sklyanin [14]), of the form 



fii^^^Hi, ...,H n ) = 0, i = l,...,n, det 



y^O, (2.19) 



which guarantee the solvability of the appropriate Hamilton-Jacobi equations and involutivity 
of Hi. A special case, when all separation relations (2.19) are affine in Hi is given by the set of 
equations 

n 

Y J <t>K\^i)Hk = ^i{\,»i), i = l,...,n, (2.20) 
fc=i 

where 4> and ip are arbitrary smooth functions of their arguments, are called the Stackel separa- 
tion conditions and the related dynamic systems are called the Stackel separable ones. 

Theorem 3 [29] Let M be a generic ujN manifold and let (Hi, H n ) be a set of n functionally 
independent Hamiltonians on N . Then, the foliation defined by (Hi, ...,H n ) is separable in DN 
coordinates and the subspace spanned by (dHi, dH n ) is invariant with respect to N* if there 
exist some functions such that 

n 

N*dHi = dijdHj, i = 1, n. (2.21) 

Hence, the distribution defined by (Hi, H n ), spanned by the Hamiltonian vector fields Xn t , is 
invariant with respect to N. 

The formula (2.21) can be written in the equivalent form 

n 

didHi = ^2a ij e dH j , i = l,...,n, (2.22) 
i=i 

which will be called a quasi-bi- Hamiltonian representation of separable dynamics. Notice that 
the projection of bi-Hamiltonian chains (2.5) onto N takes the form of (2.22). Indeed from 
(2.12) we get 

9idH^ = TT lD Wdh ( f ) W=(7T lD dh < f ) )W 

m 

= (mdh?)W-J2{z k (h?)x[ k) )w 

k=l 

m 

= (nodhgjW - (z k (h^)n dhP) \ M (2-23) 
k=i 

m 

= vrok dh^W - J2 Z k(h?)W vrok ^fV 



x» ,.. y y - ■- 

k =i 

m 

9 dH^ - Y Z k (h^)W9 dH[ k) , 



k=i 

where tt\n means a plain restriction of ir to J\f. 
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Theorem 4 [29] A n-tuple (Hi, H n ) of separable functions on ujN manifold M is Stackel 
separable iff additionally to the condition (2.21) we have 

n 

N*daij = y^a ik da k j, i,j = l,...,n. (2.24) 
k=i 

Equivalently, a n-tuple (h^\...,h^) of functions from Gel'fand-Zakharevich bi-Hamiltonian 

chains on bi-Poisson manifold M is Stackel separable iff a Poisson pencil ir^ is projectible onto 

(i) 

a symplectic leaf N of ttq and h\ functions are affine in Casimir coordinates Ci,i = l,...,m, 
i.e. 

Z k (Z l (h^)=0 (2.25) 

for all k, l,j = 1, m and i = 1, n m . 



Hence, in (2.21) 
are c-independent, so 



Z k (h? ) ) = c$ (2.26) 



N*dH?=dH% 1 -Y J 0$dH { V. (2.27) 

k=l 

In such a case the Stackel separation conditions (2.20) take the form 

m 

Y J O i j (\ i ,iJi i )h^(\ i ) = ^ i (\ i , i i i ), i = l,...,n, (2.28) 

where h^\\i) is a Casimir h^\^) evaluated at £ = A« and cr* is a Casimir multiplier. For 
the majority of known integrable systems ^(A^/iJ = ^(Aj,/zJ and cr*(Ai,//J = aj(\i,fi i ), and 
then, the separation conditions (2.28) can be presented in a compact form as n copies of the 
so-called separation curve 

m 

Y,<rj&»)h U) (0 = i>(t;,»), = (Ai,tt), i = I,-,"- (2-29) 

For the uniqness of the representation (2.29) we chose the normalization condition o~ m = 1. 
Let us define the following matrix 



/ z x (0 ••• Zx (/»("») (0 \ 



a(0 



(2.30) 



\ z ro (hW(0 ••• z m (h^)(0 J 

Lemma 5 [29] : 

(i) Determinant of a(£) coincides with the minimal polynomial of N 

deta(f) = Vdet^TT/V). (2.31) 
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(ii) Casimir multipliers <7*(Aj,/Xj) are given by the relation 

a(\ i )a i (\ i ,tii) = 0, (2.32) 

where a 1 = (a\, ...,a % m _ x , 1) T and a(Aj) is the matrix a(£) evaluated at £ = Aj. 

Up to now we sketch the passage from degenerated Poisson pencil tt\ to nondegenerated 
one 0\ : tt\ — > 0\. According to the further applications, we are interested in the inverse 
passage: 9\ — > tt\, as we would like to construct GZ bi-Hamiltonian chains from a given quasi- 
bi-Hamiltonian ones. 

Lemma 6 Let 6 be a Poisson tensor on manifold M of dim A/" = m, and let it be a second 
order contravariant tensor field on M. = M ® R ® . . . ® R o/ i/ie same ranA; as 6*, i.e. w/ai/t 

fe 

additional Casimir functions ip 1 ,...,ip k . Moreover, let the restricton of it to M (the level set of 
ip 1 = ... = ip k = 0) coincides with 6. Then, the operator it is also Poisson. 

Proof. Let us complete the set of functions <p i with some functions Xj which parametrized J\f 
submanifold, to a coordinate system (x,ip) on M. Then, the matrix of the operator tt has last 
k rows and columns equal to zero, while the m dimensional upper left block coincides with 6, 
which is Poisson. ■ 

Let 6o and 6\ be compatible Poisson tensors on M. Then, let ttq and tt\d be Poisson 
extensions of #o an d #i from J\f to M, such as in Lemma 6, with a common set of Casimirs 
(Pi, <p k . Of course, ttq and tt±d are still compatible. Now, let tt\ be another Poisson extension 
of 6\ from Af to M. with a new set of Casimirs ip 1} ip k . The question is: when tt\ is compatible 
with 7T0? 

Assume that there are vector fields Zi,i = l,...,k transversal to the submanifold M, such 
that Lz^o = 0, i = 1, k, and tt\ has the form 

k 

7Ti = 7T 1D + Y,Xi AZ U (2.33) 
i=l 

where Aj are some vector fields. 

Lemma 7 A sufficient condition for compatibility of ttq and tt\ is that Aj are Hamiltonian 
vector fields with respect to ttq . 

The Jacobi identity for a tensor tt is equivalent to the relation [7r,7r]s = 0, where [., \s is a 
Schouten bracket. Hence, the compatibility of two Poisson tensors tt and tt\ is equivalent to 
the relation [7ro,7ri]s = 0. Before we start the proof we will therefore remind some basic facts 
about Schouten bracket. 

The Schouten bracket [37], [38] in general is a bilinear mapping A q (M) x A r (M) — > A q+r ~ 1 (M) 
that with every a-vector and r-vector associates a (q + r — l)-vector [Q, R]s (Schouten bracket 
of Q and R) that is skewsymmetric in the sense that 

[R,Q] s = (-l) qr [Q,R]s (2.34) 
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and that satisfies the Leibniz identity 

[Q, R A P]s = [Q, R] S AP+ (-l) (q - 1)r R A [Q, P] s , (2.35) 

where P is any p- vector £ A p (Ai). Using the above properties of the Schouten bracket one can 
show that if X, Y are some vector fields, then [X, Y]s = [X,Y] = LxY is a usual Lie bracket 
(commutator) of vector fields and if moreover P is a bivector then 

[X AY,P] S = Y A[X,P} S -X A[Y,P} S , [X,P] S = L X P (2.36) 

Proof. From properties of the Schouten bracket 

k k 

[7T0,7Tl]s = Ko^lD + y^^t A Zj}s = y^[7T ,Xj A Zj] S 

i=l i=l 

k k 

= J2( Z i AL ^o-X i AL Zi Tr ) = J2 Z i AL x i ^o- (2.27) 
i=i i=i 

Hence, the most general condition for the compatibility is 

k 

Zi A L Xi 7T = 0, (2.38) 

i=i 

with the strong solution Lx^o = 0, i = 1, k, which takes the place when Xi are Hamiltonian 
vector fields with respect to 7To, i.e. there are some functions fi, fk such that 

Xi = ir dfi. (2.39) 

■ 

Notice that when in addition we impose the normalization conditions Zi(tpj) = 5ij,i,j = 
1, k, and commutativity of Xi then iridipj = Xj and Xj are bi-Hamiltonian vector fields 

ir dfi = mdifi = Xi, (2.40) 

which is just our case (2.12) with ip i = fi = 

Finally, let us mention the problem of explicit construction of DN coordinates (A, jj). Assume 
that j\f is parametrized by a set of coordinates {vi}f2ii n °t necessary canonical. The first part 
of the transformation (A,//) — > (v), i.e. A — > v is reconstructed from the characteristic equation 
of N 



v / det(CI — N) = det a(£) = =^ A, = \i{v), i = l,...,n. (2.41) 

The reconstruction of missing part of the transformation: fi — ► v, in general case, is far from being 
trivial and is discussed in [24]- [29]. In the following we restrict ourselves to the canonical point 
transformations (A,//) — > (q,p), where reconstruction of the remaining n relations pj = pi(\,(i) 
from the first n ones qi = %(A) is standard. 

In the next Sections we apply this theory to dynamic systems on Riemannian manifolds, 
constructing all those which can be lifted to the Gel'fand-Zakharevich bi-Hamiltonian chains. 
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3 Introduction to separable Riemannian dynamics 



3.1 Hamiltonian representation 

Let (Q,g) be a Riemann (pseudo-Riemann) manifold with covariant metric tensor g and lo- 
cal coordinates q l ,...,q n . Moreover, let G := g -1 be a contravariant metric tensor satisfying 
Y^j=\9ij^ k = The Levi-Civita connection components are defined by 

1 - r) 

T )k = g E GU ^iJ + d 39kl ~ d l9jk ), 0, = (3.1) 



The equations 



j tt + r i jk 4q* = G ik d k V(q), i=l,...,n, q t = ^ (3.2) 



describe the motion of a particle in the Riemannian space with the metric g. Eqs.(3.2) can be 
obtained by varying the Lagrangian 

£{q,<lt) = \Y.9id4-V{q) (3.3) 

and are called Euler-Lagrange equations. Obviously, for G = I equations (3.2) reduce to Newton 
equations of motion. 

One can pass in a standard way to the Hamiltonian description of dynamics, where the 
Hamiltonian function takes the form 

H ^P) = iZ^- C = \iZ G ij piPj + V(q), R := = E^ (3-4) 
i=l y * i,j=l y * j 

and equations of motion are 



{'.).■-■('- 0(1)- 

t (3-5) 
i _ dH dH 

X# denotes the Hamiltonian vector field and the whole dynamics takes place on the phase space 
M = T*G in local coordinates (q 1 , q n ,p±, ...,p n ). 

Of special importance is the geodesic motion V(q) = 0, with Euler-Lagrange equations 

qlt + T) k q^ = 0, i = l,...,n (3.6) 

and Hamiltonian representation 



Q 
P 



n 

X E , E = -J2G lJ PiPj . (3.7) 



6 dE L .. 
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3.2 Separated coordinates 

We are going to present a complete theory of a subclass of one-particle dynamics, containing 
Liouville integrable, separable, (quasi-)bi-Hamiltonian systems with n quadratic in momenta 
constants of motion. Stackel [l]-[3] gave the first characterization of the Riemann (pseudo- 
Riemann) manifold (Q, g) on which the equations of geodesic motion can be solved by separation 
of variables. He proved that if in a system of orthogonal coordinates (A, /i) there exists a non- 
singular matrix ip = (^(A*,)), called a Stackel matrix such that the geodesic Hamiltonians E r 
are of the form 

n 

1=1 

then E r are functionally independent, pairwise commute with respect to the canonical Poisson 
bracket and the Hamilton-Jacobi equation associated to E\ is separable. 

Then, Eisenhart [5]- [7] gave a coordinate- free representation for Stackel geodesic motion 
introducing a special family of Killing tensors. As known, a (l,l)-type tensor B = (or a 
(2, 0)-type tensor B = (B l i)) is called a Killing tensor with respect to g if 

{ UBGfp iPj , E} 6q = or { UBfPiPj , E} 6q = 0. 

He proved [5]-[7] that the geodesic Hamiltonians can be transformed into the Stackel form (3.8) 
if the contravariant metric tensor G = g^ 1 has (n — 1) commuting independent contravariant 
Killing tensors A r of a second order such that 

E r = lY. A r j PW> ( 3 - 9 ) 

i,j 

admitting a common system of closed eigenforms a« 

(A* - viG)ai = 0, dai = 0, i = l,...,n, (3.10) 

where v\ are eigenvalues of (1, 1) Killing tensor K r = A r g (K* = gA*). 

From now on, separated canonical coordinates will be denoted by (A, /i) and natural canonical 
coordinates by (q,p). For n degrees of freedom, let us consider n Stackel Hamiltonian functions 
in separated coordinates in the following form 

Hr = \Y^v i r G ii ^ l +V r {\) = ^ T K r G^ + V r {\), r = l,..,n, (3.11) 
i=i 

where n = (/x l5 [i n ) T ,&nd V r (X) are appropriate potentials separable in (A, /u) coordinates. For 
the integrable system (3.11) the Stackel separation conditions (2.20) take the general form 

j2<f>Hmk = \M\%l+J l (X) i = l,...,n, (3.12) 
k=l 
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where /i,7j,</>f are arbitrary smooth functions of its argument and the normalization (f>f = 
1, i = 1, n is assumed. To get the explicit form of = -Hfc(A, //) one has to solve the system 
of linear equations (3.12). The results are the following 



G l 



r+1 



detVF 



detW^ 

detW ii; 



where 



W 



i=l 

/ ^(A 1 ) 0? (A 1 ) ••• 

1 / \ tin j 2 / \ n 



= ("!) 

det W ir 
detW ' 

^(A 1 ) l\ 



(3.13) 
(3.14) 

(3.15) 



and is the (n — 1) x (n — 1) matrix obtained from W after we cancel its ith row and fcth 
column. Then the Stackel matrix ip is given by 



<P 



( ^(A 1 ) 0i (A 1 ) 

MX 1 ) JiQ?) 



\ /„(A») /n(A") 



/iCA 1 ) /TIaT) 



/n(A") TW) J 



(3.16) 



In our further considerations we restrict to the case when </>f(Aj) are monomials. Actually, 
we assume that 



#(Ai) = <^(A 4 ) = (A*)" 1 *, m n = < m ra _i < ... < mi <E N, 



(3.17) 



which is a sufficient condition for the existence of bi-Hamiltonian representation of Gel'fand- 
Zakharevich type with separation conditions and the separation curve in the form 



m 1 

5>,(A>W(A*) = /(A >f + 7 (A*), i = l,..., n , 

3=1 
m 



(3.18a) 
(3.18b) 



3.3 Integration by quadratures 



Let us consider a set of coordinates {A*,/ij}f =1 on T*Q canonical with respect to 6q. One can 
try to linearize the system (3.11) through a canonical transformation (A, ji) — > (6, a) in the form 
^ = §a~'^ = f^' where VF(A, a) is a generating function that solves the related Hamilton- 
Jacobi (HJ) equations 



E r {\ 



CL r , 



0,...,n. 



(3.19) 
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The HJ equations (3.19) are nonlinear partial differential equations that are very difficult to 
solve. In general it is a hopeless task. However, there are rare cases when one can find a 
solution of (3.19) in the separable form 

n 

W{\a) = Y J W i {\i,a) (3.20) 

i=i 

that turns the HJ equations into a set of decupled ordinary differential equations that can be 
solved by quadratures. Such (A, fi) coordinates are called separated coordinates. In the (6, a) 
coordinates the flow d/dtj associated with every Hamiltonian Hj is trivial 

— = 0, - = «%, .,j = l J ...,n (3.21) 

and the implicit form of the trajectories X(tj) is given by 

dW 

b l (X,a) = -r — = dijtj + const, i,j = l,...,n. (3.22) 

As explained in the previous section, separated coordinates are defined by n relations of the 
form 

<p i (\ i ,fi i ,H 1 ,...,H n ) = 0, i = l,...,n (3.23) 

joining each pair (A 1 , //J of conjugate coordinates and all Hamiltonians Hi, i = 1, ...,n. Fixing 
the values of the Hamiltonians Hj = const = cij one obtains an explicit factorization of the 
Liouville tori given by the equations 



ip i {\\n i ,a l ,...,a n ) = 0, i = l,...,n 

t (3-24) 

^i(A ,— -pai,---,^) = 0, « = 1, 
d\ 

i.e. a decupled system of ordinary differential equations. 

Let us now solve explicitly the Hamilton- Jacobi equations (3.19) and related dynamic equa- 
tions, with respect to the evolution parameter t r , written in (A,//) coordinates when separation 
conditions are given in the form (3.12), (3.17): 

1 .„.-, fdW^ 2 



§/(A l ) ( ^- ) + 7(A*) = ^(AT 1 + a 2 (AT 2 + - + a n = a(A* 



^ W = J2Wi(Xi,a) 



2 Jy '\d\ 



i=i 

2 



V( Al ) (^T ) + '"(A') - ^(AT 1 + «2(A')"- + - . + u„ = r,(A') 

4 
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'a(0-7(0 



4 



dm tJ V2/(6K0-7(0] 



—==d(,=t r 5 ri + consti, i = l,...,n, (3.25) 
where eqs.(3.25) are implicit solutions called the inverse Jacobi problem. 



4 Systems of Benenti type 
4.1 Separable geodesies 

Among all Stackel systems a particularly important subclass consists of these considered by 
Benenti [11], [12], [13] and constructed with the help of the so-called conformal Killing tensor. 
Let L = (L l j) be a second order mixed type tensor on Q and let L : M — > R be a function on M 

n 

defined as L = \ ^ (LG)^piPj. If 

{I, E} 6o = kE, where k = {e, E} 6o , e = Tr(L), (4.1) 

then L is called a conformal Killing tensor with the associated potential e = Tr(L). If we assume 
additionally that L has simple eigenvalues and its Nijenhuis torsion vanishes, then L is called a 
special conformal Killing tensor [39]. 

For the Riemannian manifold (Q,g,L), geodesic flow has n constants of motion of the form 

^ n 1 n 

E r = 2 E 4 J PiPj = j E (K r G) ij PiPj , r = 1, n, (4.2) 

jj=i jj=i 

where A r and iT r are Killing tensors of type (2,0) and (1,1), respectively. Moreover, as was 
shown by Benenti [11], [12], all the Killing tensors K r with a common set of eigenvectors, are 
constructed from L by 

K r+1 = Y,PkL r - k , (4.3) 

k=0 

where p r are coefficients of the characteristic polynomial of L 

det((I-L)=e + PiC- 1 + - + Pn, Pb = l- (4-4) 
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From (4.3) we immediately find that 



K r+1 -LK r = p r I, -LK n = Pn I, K 1 = I. (4.5) 

Lemma 8 From (4-3) and (4-5) it follows that appropriate Killing tensors K r are given by the 
following 'co factor' formula 

n—l 

cof(ZI-L) = Y,K n - t e, (4.6) 

i=0 

where cof(A) stands for the matrix of cof actors, so that cof(A)A = (det A)I. Notice that K\ = I, 
hence A± = G and E\ = E. 

Proof. 

n—l n—l 

^I-L)(J2 K n-iC) = Y,( K "-^ +1 - LK n-iC) 
i=0 i=0 

= —LK n + (K n - LK n -i)£ + ... 
+(K 2 - LK^e' 1 + K^ n 

= /det(£/-L). 

■ 

According to the above results, the functions E r , satisfy 

{E s ,E r } eo =0, (4.7) 

and thus constitute a system of n constants of motion in involution with respect to the Poisson 
structure 9q. So, for a given metric tensor g, the existence of a special conformal Killing tensor 
L is a sufficient condition for the geodesic flow on Q to be a Liouville integrable Hamiltonian 
system with all constants of motion quadratic in momenta. 

It turns out that with the tensor L we can (generically) associate a coordinate system 
on Q in which the geodesic flows associated with all the functions E r separate. Namely, let 
(A 1 (g), X n (q)) be n distinct, functionally independent eigenvalues of L, i.e. solutions of the 
characteristic equation det(£I — L) = 0. Solving these relations with respect to q we get the 
transformation A — > q 

^ = ai(A), i = l,...,n. (4.8) 

The remaining part of the transformation to the separation coordinates can be obtained as a 
canonical transformation reconstructed from the generating function W(p, A) = ^2iPiCti(X) in 
the following way 

Vi = dW g% X) =»Pi = A(A,M) i = l,.-.,n. (4.9) 
In the (A, p) coordinates the tensor L is diagonal 

L = diag(A 1 ,...,A n )^A n , (4.10) 
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while the geodesic Hamiltonians have the following form [22] 

E '—\t^^l r = l,....n, (4.11) 

i=l 

where 

a 4 = n ( Ai - Afc )' ( 4 - 12 ) 

fc=l,...ra, fc^i 

p r (X) are symmetric polynomials (Viete polynomials) defined by the relation 

n 

det(£I - A) = (£ - A 1 )^ - A 2 )...(e " A n ) = ^ p r r" r (4.13) 

r=0 

and fi are arbitrary smooth functions of one real argument. From (4.11) it immediately follows 
that in (A, n) variables the contravariant metric tensor G and all the Killing tensors K r are 
diagonal 

G ij = Up_ 5 ij , K j = -^L S ). (4.14) 

Remark 9 When fi{X l ) is a polynomial of order < n the metric is flat, if the order of f is 
equal n+1 the metric is of constant curvature. 



Moreover, as 



A*d\ i = X i dX i , K* r dX i = -^dX i = vi +1 dX\ 

dx l + 



then multiplying both sides of eq.(4.15) by G we get 

(GK; - v l r G)dXi = 0^{A* r - vlG)dXi = 0, i = 1, n, 



(4.15) 



(4.16) 



i.e. the tensorial Eisenhart realization (3.10) of the Stackel results. 

To ensure that in the (A,/x) coordinates the geodesic Hamiltonians (4.11) are separable it is 
sufficient to observe that in these coordinates they actually have the Stackel form (3.8) with the 
related Stackel matrix 



/ (A 1 )"' 1 (A 1 )"' 2 
' /i(A x ) h(X r ) 



V 



( A n ) n 

My 



(A") n ~ 2 

/n(A") 



(4.17) 



/n(A") /n(A") 



A n ) / 



It means that Stackel separation conditions (3.12) take the particular form 



E 1 (XT- 1 + E 2 (XT~ 2 + ... + E n = ^(A>f, 



1, ...,n. 



(4.18) 



For fi(X l ) = f(X l ) eqs. (4.18) can be represented by n different copies (£, /j,) = (A 1 , /itj, i = 1, n 
of some curve 

EiC' 1 + E 2 C~ 2 + + E n = \f(i)^ (4.19) 
called the separation curve of geodesic motion. 



17 



4.2 ( Quasi-) bi-Hamiltonian chains 

The special conformal Killing tensor L can be lifted from Q to a (1, l)-type tensor on M = T*Q 
where it takes the form 

The lifted (1, 1) tensor N is Nijehuis torsion free, like the L one, hence it defines a recursion 
operator on N . An important property of N is that when it acts on the canonical Poisson tensor 
6?o it produces another Poisson tensor 

* = Wo=(_°, (4.21) 

compatible with the canonical one (actually 9q is compatible with N k 8 for any integer k) making 
N the ojN manifold. It is now possible to show that the geodesic Hamiltonians E r satisfy on 
M = T*Q the set of relations [32] 

Q\dE r = 0odE r+ i — p r 0Q(lEi, E n+ \ = 0, r = 1, ...,n. 

t (4-22) 

N*dE r = dE r+1 - Pr dE l7 N* = 0~ 1 9 l = 



L* -F 
L 



which are a particular case of the quasi-bi-Hamiltonian chain (2.22) [40]. Obviously {E s ,E r }g 1 = 
0, as follows from (4.7) and (4.22). Notice that relation (4.5) is immediately reconstructed from 
(4.22). Indeed, according to 



r 



^,..,^J =KG P , P = ( Pl ,..,Pn) T , (4.23) 
dp components of (4.22) are 

LK r Gp = K r+ \Gp — p r Gp 

(4.24) 

LK r = K r+ \ — p r I. 

On the other hand 

LK r Gp = K r+ \Gp — p r Gp 

4- (4.25) 
L r = ^p T LK r Gp = ^p T K r+1 Gp - ^p r p T Gp = E r+1 - p r E u 
so, for r = 1 we get the relation 

L = L 1 =E 2 - p x E x (4.26) 



18 



and immediate verification that L is a special conformal Killing tensor, as 

{L,E 1 } ea = {E 1 , Pl } 0Q E 1 = kE u Pl = —TrL. (4.27) 

As was presented in this section, tensors L* and A* have a common set of closed eigenforms 
(4.15), (4.16), hence obviously 

(G (s) iC - vlG {s) )d\i = (4 s) * - v l r G {s) )dXi = 0, i = 1, n, (4.28) 

where 

4 s) = L s A r = K r L s G = K r G {s \ s = ±1, ±2, ... . (4.29) 
It means that are contravariant Killing tensors of the metric G^ and functions 

£7<>) = ^p T A^p, r = 1, n (4.30) 

are in involution. Moreover, L is the special conformal Killing tensor for each G^ s \ Indeed, from 
(4.24) we have 

L s+1 K r G =L s K r+1 G - p r L s G 
LK r L s G = K r+l L s G - p r L s G 
LK r G {s) = K r+1 G {s) - Pr G {s) 

4 S) =E { r i 1 -p r E{ S \ (4.31) 



where 



4 s) = ^p T LA^p. (4.31a) 



Hence, for L = L\ , we find the condition (4.27) 



{L {s \e{ s) }0 {) = {E^.p^E^ = K E{ S) . (4.32) 



Let us denote by G^ a basic flat contravariant metric, which in A coordinates takes the 
form 

(gWJ J = ±-5 ij , (4.33) 



i.e. fi(X l ) = 1, i = 1, ...,n (4.14). It means that in the appropriate separation curve for geodesic 
motion (4.19) /(£) = 1. Moreover, the metric tensor G for which the separation curve (4.19) 
contains /(£) in the form of Laurent polynomial 

i=m,2 

f(0=^2 a iC, mi,m 2 eN, cii = const, (4.34) 
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is constructed from the basic metric in the following way 

G = /(L)G (0) . (4.35) 

Now we pass to the bi-Hamiltonian representation. On the extended phase space M = 
T*Q x R, the extended geodesic Hamiltonians 

e r = E r + cp r , r = l,...,n, eo = c, (4.36) 

where c is an additional coordinate, satisfy the following bi-Hamiltonian chain [32] 

iTodeo = 



(4.37) 



with the Poisson operators ttq and tti 



#0 












7T1 







-K\de\ 




nide n -i 




iride n 




0i 


Oodei 


-(Oode 1 f 






(4.38) 



Both Poisson tensors 7To and tt± are compatible and degenerated. The Casimir of ttq is eo and 
the Casimir of tt\ is e n and they start and terminate the bi-Hamiltonian chain (4.37). Of course 
{e s ,e r } no = {e s ,e r } ni = 0. The projections of 7ro,7ri onto a symplectic leaf of ttq (c = const) 
along Z = J|; reconstructs our nondegenerate Poisson tensors 9q,9\. If we introduce the Poisson 
pencil 7Tg = -K i — ^ttq, the chain (4.37) can be written in a compact form 

n 

7r c <fe(£) = 0, e(0 = ^e„_ r e r , (4.39) 

r=0 

where eg is a Casimir of the Poisson pencil 7rg, depending polynomial on £, i.e. we deal with a 
Gel'fand-Zakharevich type system. The projection of (4.37) onto symplectic leaf of ttq : c = 
reconstructs the quasi-bi-Hamiltonian chain (4.22) on A/" in a compact form 

9td(E(S))+p(S)9 d(E 1 ) = 

t (4-40) 
(N*-ZI)d(E(O)+p(Od(E 1 ) = 0, 



where 



n-l 



E(0 = En-rC, P(0 = det(£/ -L) = J2 Pn-MC- (4.40a) 

r=0 r=0 

If we start with a bi-Hamiltonian chain (4.37) or a quasi-bi-Hamiltonian chain (4.22) written 
in the 'physical' (original) coordinates (q,p) then we can usually find the functions q(X) that 
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constitute the first half of the transformation (X, p) — > (<z,p) simply by taking the functions 
p r (q) from the Hamiltonians (4.36) or directly from (4.4) and solve the system of equations 
p r (q) = p r (A) , r = 1, . . . , n with respect to q. The remaining functions p(X, p) are calculated 
according to (4.9). 

4.3 Separable potentials 

What potentials can be added to geodesic Hamiltonians E r without destroying their separability 
within the above schema? It turns out that there exists a sequence of basic separable potentials 

(k) 

V; k = ±1,±2,..., which can be added to geodesic Hamiltonians E r such that the new 
Hamiltonians 

H r (q,p)=E r (q,p)+V r ^(q), r=l,..,n, (4.41) 

are in involution with respect to 9q and 9\ and are still separable in the same coordinates (A, p). 
It means that H r follow the quasi-bi-Hamiltonian chain (4.22) 

N*dH r = dH r+1 - p r dH 1: H n+1 = 0, r = 1, n, (4.42) 

while for potentials we have 

L*dV r = dV r+ i -PrdVu r = l,...,n, (4.43) 

t 

r 

dV r+1 = Y,Pk{L*) r - k dV 1 

k=0 

t (4-3) 

dV r+1 = K; +1 dV! => d(K; +1 dVi) = 0. (4.44) 
The relations (4.44) were derived for the first time by Benenti [11], [12]. 

Theorem 10 The basic separable potentials Vr™^ are given by the following recursion relation 
[23], [27] 

y r (-+D = V^l-p r vi m \ (4.45) 

and its inverse 

Pn 

Proof. The proof is inductive. We show it for positive potentials. Assuming that potentials 
V r {m) fulfil condition (4.43), we prove that potentials V^ m+1 ^ fulfil the same condition. The 
condition (4.43) is true for the first nontrivial potentials = p r , which are coefficients of 
characteristic polynomials of a special conformal Killing tensor L [32] 

dpr+i = L*dp r + p r dpi, r = l,...,n. (4.47) 
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Notice that (4.47) is a special case of the condition (2.24) and follows from the Stiickel separa- 
bility. Then we have 

L * dV (m+l) + p rdV f™+l) 

= L*d{V™ - p r v} m) ) ~ Prd{vt ] - PiV} m) ) 

= L*dV^l - p r L*dv} m) - V} m) (L*dp r + Pr d Pl ) 

+ Pr {dvi m) ~ Pidvi m) ) 

= L*dV^l - Pr L*dv} m) - V} m) dp r+1 + Pr L*dv} m) 

= L*dV%l - Vi m) dp r+l = dV^l - Pr+ M m) ~ V} m) dp r+l 

= d{V<$> - Pr+ M m) ) = dv}™ +1 \ 

m 

Notice that = —5 rn - m ,m = 0, ...,n — 1, Vr^ = p T ^Vr^ 1 ' ) = Such a notation will 

Pn 

be useful in the case of deformed Benenti systems. 

Now, again the extended Hamiltonians li r :JVxl^l 

h r = H r + cp r (4.48) 

satisfy the bi-Hamiltonian chain (compare with (4.37)) 

nodho = 

-Kodhi = X\ = iridho 

TTodIl2 = X2 = 7T\dhi 



(4.49) 



TTodh n = X n = TTidh n -i 
= lTidh n 



with ttq as in (4.38), 



^ -(M«>i) T J 

and {h s , hr}^ = {h s , h r } ni = 0. If we use the following notation 

n-l 
3=0 

then the recursion formulas (4.45) and (4.46) can be written in the compact form [41] 

V(0 (k+1) = XV{0 (k) - det(£/ - L)v} k) (4.51) 

and 

K (o<-*-i> = I (v(0^ - det Z~L L) v ^ k) ) ^ 
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and our bi-Hamiltonian chain (4.49) is given by 

ir ( dh(O = 0, h(0 = H(0+cp(0, (4.53) 
while the corresponding quasi-bi-Hamiltonian (4.42) chain takes the form 

9^dH(0 + piWodHifap) = 0, p(0 = ^ 

t (4-54) 

(N*-ZI)dH(Z)+p(Z)dH 1 (q,p)=0. 

In (A, p) coordinates the full (i.e. with a non-zero potential part) Hamiltonians (4.48) of bi- 
Hamiltonian chain (4.49) attain the form 

h r (X,p,c) = -^j + cp r (X), r = l,...,n. (4.55) 

i=i 

Lemma 11 Nontrivial potentials V^ n 1+k ^ and Vr ^ k = 1,2, ... enter the separation curve 

H^ 1 - 1 + H 2 C~ 2 + ... + H n = l -f(i)p 2 + 7 (0 (4.56) 



as 7 (f) = -C~ 1+k , ~C k , hence 7i (A*) = -(A*) n ~ 1+fc and 7i (A*) = -(V)" fe , respectively. 

roof. Potentij 
Killing tensor L 



Proof. Potentials = p r are coefficients of characteristic equation of the special conformal 



Then, we define y( n+fc ) potentials by a generating equation 

C n+fc + V^+^f- 1 + ... + y n (n+A;) = 0. (4.58) 
The recursion formula (4.45) is reconstructed as follows. From (4.58) we have 

Elimination of £ n via (4.57) leads to the form 

c+k+l + (v , ( „ +fc) _ + + (v K„+ fc ) _ p^^/^))^ _ Pn v{ n+k) = 

4 (4.60) 

T/ (n+fe+l) _ T /(n+fc) T /( n+fc ) 
v r — v r+1 p r .i/ 1 

For the inverse potentials we have from (4.57) 

J_^-i + ^-2 + ... + ^tzl + r i =0 , (4.61) 

Pra Pn Pn 
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so 

y( 1) = rj,^ 

we define fc ) potentials by a generating equation 

Fi (-*)^n-l + ... + y(-k) + £-* = Q _ (4>62) 

The recursion formula (4.46) is reconstructed as follows. From (4.62) we have 

v (-k)^n-2 + + Fn (-*)^-i + = o. (4.63) 

Substituting (4.61) to eliminate £ _1 we get 

y{-k) e -2 + + y(-k) { _]_ e -l _ _ Pn z l ) + r fc-l = q 

Pn Pn 

^ (4.64) 

" r " r— 1 n 

rn 

■ 

5 One-hole deformation of Benenti systems 

Separable systems on Riemannian manifolds considered by Benenti belong to important but 
very particular subclass of such systems. In this context, a question about classification of 
all separable systems on Riemannian manifolds, with n quadratic in momenta constants of 
motion, arises. The classification can be made with respect to the admissible forms of Stackel 
separability conditions. The right hand side of the conditions (3.12) is always the same for the 
class of systems considered 

r.h.s. = l -h{\ l )tf + 7i (V) = Hi), (5.1) 

so different classes of separable systems are described by different forms of the l.h.s. of Stackel 
conditions. In the simplest Benenti case, it is given by the following polynomial form 

l.h.s. = H 1 C~ 1 + H 2 e- 2 + ... + H n , £ = A 1 ,...,A n . (5.2) 

We will show, that all other classes, given by some polynomial in A, are appropriate deformations 
of the Benenti class. 

First, let us define a 1-hole deformation of the Benenti class. Consider the following separa- 
bility condition 

H^n+l)-! + # 2 £(n+l)-2 + + jj n+1 = ,/,(£, /i), H ni = (5.3) 

and the Benenti separability condition with the same ip representation 

HiC' 1 + H 2 C- 2 + ... + H n = n). (5.4) 
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Notice that all Benenti systems are classified by different forms of ip, i.e. by fi(X l ) and j^X 1 ). 
A missing monomial (a hole) in (5.3) is H ni ^ n+1 ^ ni . Using the characteristic equation of a 
conformal Killing tensor L of the Benenti system (5.4) 

e+p 1 e- 1 + ...+p n = o, 

for the elimination of £ n , equation (5.3) can be transformed to the form 

(H 2 - p^)^ 1 + + (H n+1 - Pn H x ) = ^,P) (5-5) 
and hence, comparing it with (5.4) we find 

H r = H r+1 - p r Hi, r = l,...,n, (5.6) 

with the inverse 

H r+1 = H r -^—H ni ^, r = 0,...,n, (5.7) 
Pm-l 

where Ho = 0, p = 1. 

Notice, that formula (5.7) applies separately to the geodesic and the potential parts, i.e. 

E r+ \ = E r — — E ni -i, (5.8a) 

Pni-l 

Vr+i = V r --^V ni - 1 , r = 0,...,n. (5.8b) 

Pni-l 

5.1 Geodesic part 

Let us first look onto geodesic Hamiltonians 



E r = (p m -l K r-l - Pr-l K ni-l)—^—Gp, T = 1, ...,n+ 1. (5.9) 

z Pm-l 



Using the known relation for the Benenti chain 

p r I = K r+ i - LK r 

we get 

1 TP . - 1 1JT 



Ei = --^-£ ni -i = -j^-Ar K ni -iG V 

l J n i — l In i — l 

= \P T G P ==>■ G = ——^—K ni -iG 



(5.10) 



p ni -i 



and _ 

E r = \P T \j^rS K n,K r -l ~ K ni ^K r )G}p 



(5.11) 



= \p T K r G p^K r = K r - K r ^K ni (K ni ^)- 1 . 

Although we know from the construction that E r are in involution, as they are Stackel geodesies, 
but here we prove it in a coordinate free way. 
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Lemma 12 K r are Killing tensors of the metric G which share the common set of eigenfunc- 
tions, i.e. E r are in involution. 

Proof. From involutivity of Benenti Hamiltonians we have 

= {H r ,H s }e = {E r + V r ,E s + V s } eQ = {E r ,V s }g + {V r , E s } e „ 



(5.12) 



{Er, V s }e = {E s , V r } do => {E r ,p s } dQ = {E s ,p r } eo . 

Then, 

{E r+1 ,E s+1 }g = {E r , ^—E ni ^i} e() + { -?—E ni -i,E s }o 

Pm-i Pni-l 

+ { ^ E ni -i, — E ni -i}g Q 
Pn\—\ Pni—1 

= { E r,P ni -l}e 2 S E rn-1 + {Pni-1> E s}o ~ Eni-1 
Pni-1 Pni-1 

+{p r i E n 1 -l}o 2 S E ni-1 + {E ni -1, P s }o ~ 
Pni-1 Pni-1 

= 0. 



5.2 Deformed potentials 

Let us analyze the basic deformed potentials. The first potentials are the following. From (5.8b) 
and the Benenti potentials we have Vr™^ = —5 r -i jn - m for m < n + 1, m ^ (n + 1) — m and 

pr(n+i)-m = A=L, v;(-+i) =/9r ._^zl^l,... . (5.13) 
Pni-1 Pni-1 

Notice that = for m > n + 1 and V n ( " +1)_ni = 1. 

Lemma 13 Nontrivial basic potentials V^ n+1 "-^ y r ( n + fc ) an( ^ k ) ^ = i ; 2, ... enier i/je sep- 
aration curwe 

^ir + ^2r _1 + - + ^n+i = ^/(6p 2 + 7(0, fl"m=0, (5.14) 

0S7 ^) = _^(n+l)-ni ) _^+fc j _ r fc. 

Proof. We will show the following generating equations for the potentials considered 

C +k + Vi n+k \ n + .. + Of* 1 ""! + ... + V^X k) = 0, A; = 1,2,..., (5.15a) 

r fe + ^i (_fc) f + ... + 0e n+1 " ni + ... + V r [~ k) = 0, fc = 1, 2, .., (5.15b) 

y(n+l- ni ) c + + e „ + l_ ni + + V^n+l-m) = Q ( 5 .1 5c ) 
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r + v^H-piC- 1 - - - pj + - + ri+l = o 



r + (vt ] - pM^' 1 + - + (v^l - Pn v{ m) ) = o 



For the first two equations we have (m > n or m < 0) 

r+^ (m) r+-+^ ( 3 = o, K ( r ) = o 

T/ i+i 

+ + + y{rn) = q 

which reveal the known deformation relations (5.8b) 

= ^ _ Pr v} m) ^ V^l = - -^-V^\ (5.16) 

Pni-l 

between nontrivial basic potentials from Benenti class and respective deformed potentials. For 
the last case (5.15c) we have 

t 

y(n+l-m) 1 ^-(n+l-ni) 

r + — rC" 1 + ... + ^7-; — -r +1 ~ ni + ... + — r = 

y(n+l-ni)^ y(n+l-ni)^ y(n+l— m) 

t 

e+PiC- 1 + -+Pn = o 

y(n+l)-m _ Pr-l 
Pm-i 

which is a special case of deformations (5.16) related to the trivial Benenti potential y( n+1 ~ ni ) = 

0~r,n\ — 1- ^ 

Alternatively, the basic potentials W m ) , m > n + 1 can be constructed recursively as in the 
Benenti case. 

Lemma 14 The basic separable potentials , m > n + 1 are given fry the following recursion 
relation 

y r (™+D = v£l - V r ^v} m) - V r ( n+1 - n ^% (5.21) 
where V^ n+1 ^ ni and V^ n+1 ^ are given by (5.13). 
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Proof. The potentials Vr n+1 ^ ni enter the separation curve in the form (5.15c), while the 
potentials Vr a+l \ Vr m \ V^ m+1 ^ enter the separation curve in the form (5.15a) with k = l,m — 
n, m + 1 — n. The recursion formula (5.21) is reconstructed as follows. Multiplying equation 
(5.15a) for k = m — n by £ we have 

Substituting £ n+1 from (5.15a) for k = 1 and ^ n+1 ^ n i from (5.15c) we get 

+ ^^(n+l-nO^n _ ... _ ^-0) + ... + ^ = Q 

A comparison with the separation curve for the potential V^ m+1 ^ (eq. (5.15a) with k = m+l — n) 
reveals the formula (5.21). ■ 

Of course, all Hamiltonian functions H r are in involution with respect to Oq as the proof of 
Lemma 6 is valid when we change E{ — > Hi and Ei — > Hi. 

5.3 Quasi-bi-Hamiltonian representation 

The invariance of the subspace spanned by (dH\, dH n+ i) with respect to N* (2.21) is proved 
in the following theorem. 

Theorem 15 Hamiltonian functions H r fulfill the following quasi-bi- Hamiltonian chain 

dH r+1 = N*dH r + a r dHi + /3 r dH ni+1 

t (5-17) 
6()dH r +\ = 9\dH r + ot r Q§dH\ + P r 9odH ril -^i 1 
where a r = V r {n+1 \(3 r = V^ +1) ~ ni . 

Proof. We use the property of the Benenti chain 

dH r+1 = N*dH r + VWdHx, (5.18a) 

d Pr+1 = L*d Pr + V r {n) d Pl , V} n) = p r , (5.18b) 
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and recursion relations (5.7), (5.8b). Hence, we have 
r.h.s.(h.n) = N*dH r + a r dH x + (3 r dH ni+l 

= N*d{H r _, - -^-ir ni _!) + ( Pr - ErzlEm.) d ( — l — Hi 

Pm-l \ Pm-l / Pm-l 



i Pr-l (tt Pni u \ 

H \ H m -Hni-lJ 

Pni-l Pm-l 



Pm-l Pni-l 

H ni -iN*dp r _ 1 + ^ r 1 di^ni — ^ r 1 H ni -idp m 

Pni-l Pm-l Pm-l 

H ni -l 



p r d 



Pm-l 



= JV*dff r _i + p r _ l0 Lffi - / 9„ 1 _ 1 F ni _ 1 iV*d(^^) 

Pni-l 

— H ni ^ 1 N*dp r _ 1 - ■^^-H ni _ 1 dp m - p r d ( 

Pni-l Pm-l V Pni-l 

V Pm-l / Pm-l 
' -H^N^dp^ 



Pm-l 



dH r - p r d ( " ftl ) - ^ 1 H ni ^ 1 dp l — H ni ^ 1 N*dp r _ 1 

Pm-l 



dH r ~ P r d f ) - lUiiZi (TV* cZp^ + p r _!d Pl ) 





-1 


Pni- 


-1 


#ni 


-1 


Pni- 


-1 




-1 







Pm- 


1 


^m 


-1 


Pni- 


-1 


#n r 


-1 


Pni- 


-1 



= — p r d I I - — dp r . 

\ Pm-l J 

= d[H r — i? m -i J = dHr+i = /As. (5.17) 

\ Pm-l J 

m 

Notice that from (5.17) it follows that also {H r , H s }q 1 = is valid. Of course, formula (5.17) 
works separately for E r and V r in the form 

dE r+1 = N*dE r + a r dEi + (3 r dE ni+1 , (5.19a) 

dVr+i = L*dV r + a r dVi + f3 r dV ni+1 . (5.19b) 
Repeating recursive transformation (5.19b) we end up with the equivalent form 

dVr+i = A r dVi + B r dV ni+1 , (5.20) 

where 

A r = k; +1 - -^k;, B r = —^—k;. (5.20a) 

Pm-l Pni-l 
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In analogy to the Benenti case (4.23), (4.24), the d p part of (5.19a) gives us immediately the 
analog of formula (4.5) for the one-hole case 



K r+1 = LK r + a r I + (3 r K ni+ i. 
Let us introduce L function as it was done for the Benenti case 



Then, from (5.21) we find 



and 



L = - P I LGp. 



L = E 2 - a x E x - ^E ni+1 



(5.21) 

(5.22) 

(5.23) 
(5.24) 



{L,E 1 } 9{> = {E 1 ,a 1 } 9a E 1 + {E ± , /^K^+i = k°E 1 + ^E ni+l . 
Thus, for G, L is not a conformal Killing tensor. 

5.4 Bi-Hamiltonian representation 

Adding two Casimir coordinates ci,C2 (with respect to 9q) and extending Hamiltonians H r to 
the form 

h r = H r + a r c\ + (3 r C2 (5.25) 

one can put the quasi-bi-Hamiltonian chain (5.17) from T*Q into a bi-Hamiltonian chain on 
T*Q x R x R, being a composition of two bi-Hamiltonian sub-chains 

7Todh = 

TTodhi = X\ = nidho 



7Todh ni -l = X ni -i = TTldh ni -2 

= 7?id/i ni _i 



(5.26a) 



TT dh ni = 

TTO^ni+l = -Xni+1 = ^\dh ni 
7Todh n+ i = X n+ i = TTidh n 

= Tiidh n+ i , 



where /iq = c±, h ni = c 2 and 



#0 


\ 









J 












9odhi 9odh m +i ^ 


7Tl = 


-{9 Q dh{) T 
V -(9 dh ni+1 ) T 






(5.26b) 
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7To and 7Ti are compatible Poisson tensors on T*Q x E x E as in (q,p,c) coordinates tt\ tensor 
has the form (2.33) 

k 



i=i 



where 



kid 



9i 


\ 


















, Zi = (0,...,0,l,0) r , Z 2 = (0,...,0,0,1) T , 

X 1 = TTodh,! = (9 dhi,Q) T , X 2 = TT dh ni+1 = (6 dh ni+1 , 0) T 

and vector fields Zi and Xi fulfil condition from Lemma 7. The reduction of both chains onto 
the symplectic leaf c\ = c 2 = of tt q along Z\ = , Z 2 = reconstructs immediately the 
quasi-bi-Hamiltonian chain (5.17). 

If we introduce the Poisson pencil tt^ = tt\ — £ttq both chains (5.26a), (5.26b) can be written 
in a compact form 

^dh(0 = 0, h(0 = (hW(0,hW(£)f, (5.27) 



where 



h {1 \o = hoc 1 - 1 + ... + v-i, h^(o = h ni e +1 - ni + ... + k 



+i 



(5.27a) 



and the operation is applied to each component of the vector. The quasi-bi-Hamiltonian chain 
(5.17) takes the form 

p(0 + 5(o(w%) =0 



t 



(N*-ZI)dH(O+a(Z)dH (1) =0, 



(5.28) 



where 



— #i — £#o, #(i) — (Hi,H ni+ i) T , 



H (1 »(0 = fftf" 1-2 + + »»,-!, »?> = ff„, + i{"""' + ... + 



and 



5(0 



9ci dci 



\ 9c 2 dc 2 I 



(5.28a) 



(5.28b) 



The Stackel conditions for extended Hamiltonians h r take the form (3.18b) 

Cl ^(-+l) +^^(n+l)-l + ... + ^ m _^(n+l)-(n 1 -l) 
+C2 £(»+l)-"i + ^ ni+1 ^("+l)-(«H-l) + ... +^ n+1 



r +2_ni /» (1) (0+^ (2) (0 = ^/(0^ 2 +7(e), 



(5.29) 



where (£, (i) = fa), i = l,...,n. 
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6 /c-hole deformations of Benenti systems 
6.1 Deformation procedure 

Here we extend the results of the previous section onto the general fc-hole case. Let us start 
with the separability condition 



Hl ^n+k)-l + H ^(n + k)-2 + + ^ = ^ ^ 



(6.1) 



with fc-holes in £(n+k)-m ^(n+k)-n 2 ^ _ ^(n+k)-n k ^ 1 < m < ... < n k < n + k, k £ N, i.e. 
H ni = H ri2 = ... = H„ k = 0, and the separability condition for Benenti systems with the same 

H^- 1 + H 2 i n - 2 + ... + H n = il>(£, fi). (6.2) 

As for the basic potentials 

substituting this relation to (6.1) for ^( n+fe )~ 1 ) ...,£ ra we get a deformation of the chain (6.1) to 
the Benenti case (6.2) 



H r = H r+k - V^-^H! - V^ n+k - 2) H 2 



(6.3) 



where H m = ... = H n = and V?™^ are appropriate basic Benenti potentials. 



Lemma 16 Deformation of the Benenti case (6.2) to the chain (6.1), i.e. the inverse formula 
to the (6.3) one, is given by a following determinant form 



Hr-k Pr-l 
Hni-k Pni-1 



Pr-k 
Pnj—k 



Hn k -k Pn k -1 ' " " Pn k -k 



Pm-i 



Pnx-k 



(6.4) 



Pn k — \ ' ' ' Pn k -k 

Proof. First, we select from (6.3) k equations containing H ni , ...,H Hk 

tt _ jAn+k-1) 7t t.M) 7t 

H ni _ k - -V k Hi-...-V k H k , 



H„ 



_ v (n+k-l)^ 



T/W TT 

V n k -k H k- 



The solution with respect to H%,% = 1, k is given by a determinant form 

S- Wi . 

Hi = —, i = l,...,n, 
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where 



and 



W = {-if 



y(n+k 

m—k 


- 1 ) 




y(n+k 
n h ~k 


-1) 


^n k -k 



Wi = (-1) 



k+i 



Hn k -k 



y(n+k- 
n\—k 


1) 




y{n+k- 
n k ~k 


1) 


■ v (n \ 

n k -k 



with the column (V r 



(n+k—i) 

m—k ' 



Ki^-fc m i ssm S- Substituting this result to (6.3) we get 



H r „ k W + V^t k X Vi + ... + V^ k W k 

w 



H r - k V 



H ni - 



(n+k-l) 
-k 

(n+k-l) 



V 



(n) 



k v • 

K m—k 



V 



-k 
(n) 



m—k 



y{n+k-l) 
n k - u 



V 



(n) 



n k - 



y(n+k- 
n\—k 


1) 


• vitk 


y(n+k- 
n k -k 


1) 


^n k -k 



H r -k Pr-l 

H ni -k Pm-i 

H n k -k Pn k -1 



Pr—k 
Pn\—k 

Pn k —k 



Pni-1 



Pn\—k 



Pn k — 1 ' ' ' Pn k —k 
(n) 

The last step is valid due to the fact that V { = Pi, the form of the recursion formula for Benenti 
basic potentials (4.45) and the properties of determinants. It allow us to replace the arbitrary 
potential v(: a + k -' 1 ) i n determinants by the = p one. For each recursive step we have 

r (n+k— i) 



V, 



m- 



Pni—k 



(n+k- 



n k -k 



V 



Pn k —k 
(n+k—i 



-1) 



m—k+l 



Pm-kYx 



(n+k-i-1) 



^An+k— i— 1) -r An+k— i— 1) 

V n k -k+l ~ Pn k -k v l 



Pn\—k 



Pn k —k 
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V 



(n+k-i-1) 



n\— fc+1 



Pn\—k 



r (n+k—i—l) 
rik- 



v -< t -h+l Pn k -k 



The formula (6.4) applies separately to the geodesic and the potential parts. 
6.2 Deformed geodesic motion 

Let us first look onto n geodesic Hamiltonians E r , r = 1, ...,n + k, r / m, ...,nfc. Introducing 
the abbreviation 

Pni-l ' 



p(m, ...,n fc ) = 



one finds 



f 



E r -k Pr-l 
E ni -k Pni-1 



E n k -k Pn k -1 ' " " 

Using the known relations for Benenti chain 

p r I = K r+ \ — LK r 
and the property of determinants we get 



E r 



Pn\—k 

Pn k —k 

Pr-k 
Pn\—k 

Pn k —k 





Kr-k 


Pr-l ■ 


Pr-k 


1 rp 

2 P 


K-ni—k 


Pm-l ' 


Pn\—k 




K nk -k 


Pn k -1 ■ 


Pn k —k 




K r -k 


K r ■■ 


Kr-k+1 


1 T 

2 P 


K-ni—k 


K ni •• 


Kni-k+1 




K nk -k 


Kn k ■ 





-Gp 



-Gp 



(6.5) 



(6.6) 



(6.7) 



K 



ni 



Kr-k+l K r -k 



K, 



n k —k+l E-n k —k 



r 



-Gp 



= (-l) k lp T (K r D - K^D, + ... + (-l) k K r . k D k )-Gp, 

Z if 
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where 



Di = 



K„ 



K„ 



K, 



m—i+i K ni —i—\ 



K, 



n\— k+1 



*-n k ' ' ' ls -ri k — i+l Kn k —i—l ■ ■ ■ i\ nfc _fc^i 

and K m in determinant calculations are treated as symbols not matrices. 
Then, 



, i = 0, k 



Ei = 1 -P T G P 



G=(-l) k -D G, 



and 



where 



E r = -p T K r Gp ^K r = K r - Kr-iDiDy 1 + ... + {-l) k K r _ k D k D \ 



K, 



n\ — \ 



K, 



n\—k 



Kn k -l ' ' ' K-n k -k 



(6.8) 
(6.9) 

(6.10) 



Again we know from the construction that E r are in involution, as they are Stackel geodesies. 
Here we present the proof in the coordinate free form but on the level of full Hamiltonians 
(geodesic motion is included as a special case). Let us first consider basic potentials. 



6.3 Basic deformed potentials 

From (6.4) the deformed potentials are 

y(™) 



<p(m,...,n k ) 



V r-k Pr-l 

V (m) n 
v m~k Pm-i 



r (m) 



V ~;-k Pn k -1 



Pr-k 
Pri\—k 

Pn k —k 



(6.11) 



so using the recursion formula for the Benenti potentials and the properties of the determinants 
we have = —5 r _k,n-m, m < n + k, m ^ (n + k) — rii, i = 1, k, 



y^n+k)-rii _ 



Pr-l 
Pni-1 ' 



Pr—k 
Pni—k 



with the row (p„._i, -,p ni -k) missing, 



y(n+k) = ±_ 



Pr Pr-l 
Pn\ Pn\ — 1 



Pn k —1 ' ' ' Pn k —k 



Pr-k 
Pn\—k 



(6.12) 



J n k 



Pn k — 1 ' ' ' Pn k —k 



(6.13) 



Notice that Vn™^ = for arbitrary m > n + k and V^™ + ^ ™ 3 



"iji h 3 — 1) ■•■) k- 
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As in the 1-hole case, one can shaw that nontrivial basic potentials Vr n+k n% \i = 1, k, Vr n+k 



and Vr ^ / = 1,2, ... fulfill the following generating equations 

e+ k-i+i + ^(n+fc-i+i)^(n+fc)-i + ... + v^_ + k k - 1+l) = 0, (6.14a) 

r' + Vt l) t {n+k) - 1 + - + V; ( ;? = 0, (6.15b) 

^n+k-r H + yfa+fc-rai)£(n+k)-l + ^ + ^(n+fc-n;) _ ^ (6.14c) 

and hence enter separation curve 

as 7 (0 = -$("+*)-»*, -£ n+k , -C h , where i = 1, k, 1 = 1, 2, ... . 

Also as in the 1-hole case, it is not difficult to shaw that basic potentials Vr™ 1 ^ ,m > n + k 
can be constructed recursively by the following recursion relation 

k 

y r (™+D = - V r ^v} m) - K {n+fc ~ n ^iSi, (6-16) 

i=i 

where v} n+k ~ ni \i = 1, k and V r {n+k) are given by (6.12) and (6.13). 



6.4 Deformed Hamiltonian functions 

For further applications we present now a recursion formula for deformed Hamiltonian functions 
understood in the following sense. Having the Hamiltonian functions for the fc-hole case we 
construct the Hamiltonian functions for (k + l)-hole case, respectively. Let us consider the 
separability condition 

^(n+fc+l)-! + -JJ^n+k+D-2 + + Jj^^ = ^ ^ ^ 
With (k+ l)-holes in ^(n+k+l)-(n l+ l) ^ ..^ ^ n +k+l)-(n k +l) ^ ^(n+fc+l)-n fc+1) 1 < ni < ... < < 

n + k + 1, fceN, i.e. H ni+ \ = ... = H nk+ \ = H nk+1 = 0, and the separability condition 

H^+k)-i + H 2 ^+ k )^ + ... + H n+k = fj), (6.18) 

with fc-holes in ^+k)- ni ^(n+k)-n 2 ^ .. ^(n+fc)-n fc) 1 < nj < ... < nj, < n + I;, G M, i.e. 
// rai = .£f n2 = ... = ff nfe = 0, and the same ip. It means that the first k holes are fixed. From 
(6.14a) with I = 1 we have 

e+ k + y(n + k)^ n+kyi + _ _ + ^(n+fe) = Q (g ^ 

Substituting (6.19) to the first term of (6.17) we get 

+(H i+1 - F. (n+fe) ^ 1 )e (n+fe) - i + ... + (^n +fc+ i - K!:t fc) ^i) = M) (6.20) 
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and hence 



where H ni+l 



H r = H r+1 -V r (n+k) H 1 , r = l,...,n + k, 



(6.21) 



Hn k +1 — H nk+1 — 0, Vni + ^ 



V^ +k) = 0. The inverse of (6.21) 



V, 



(n+k) 



H r+1 - H r ~( n+k ) H n k+1 -1 
V n k+1 -1 



(6.22) 



is our final recursion formula. The formula is applicable separately for geodesic Hamiltonian 
functions E r and potentials V r . 

Lemma 17 Functions H r are in involution with respect to the Poisson tensor 9q. 



Proof. The proof is inductive. The 1-hole case was proved in Section 3. From the involutivity 
of A;-holes Hamiltonians H r we prove the involutivity of H r Hamiltonians from (k + l)-holes 
case. From involutivity of H r Hamiltonians we have 

= {H r , H s }g = {% + V r , E s + V s } 0() = {E r , V s }e + {V r , E s }g Q 

{H r , V s }g = {H s , V r } 6o {H r , Vj n+ % () = {H s , V r {n+k) W 
Then, from (6.22) we have 



~ y( n+k ) _ 

{H r+1 ,H s+ i} 6o = < H r , - ~ (n+fc) ^wfc+i-i 



+ { 

n k +i — 1 



V r (n+k) ~ 
~y(n+k) H n k+ i-U H s 
V n k+1 -1 



y(n+k) _ y{n+k) 

~(n+k) H n k+ i-li ~{ n +k) Hn k+i-l 
k "fe+i-1 V n k+1 -1 

f^(n+fc) 

— iff T/( n + fc ) 1 Vs 
~ ^ V n k+1 -M—^ 

\ V n k+1 -l 



2 Hfik+i-l 



2 H ' 



y(n+k) ^ 

+ {^+1^1 ' H "ho ~ (n+k) H n k+1 -1 
"fe+i-1 

~ ~ T/( n + fe ) 

\ n k+1 -l J 

~ ~ y( n+k ) „ 

{ H n k+1 -i, v s n+k) )e ~ {n+k) H nk+1 -! 

V n k+ i-l 

= 0. 



T n k+1 -l + 
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6.5 Quasi-bi-Hamiltonian representation 

Now we demonstrate that the subspace span by {dH\, ...,dH n+ k) is invariant with respect to 
N* (Theorem 3). 

Theorem 18 Hamiltonian functions H r belong to the following quasi-bi- Hamiltonian chain 

k 

dH r+1 = N*dH r + a°diii + ^ a^dH ni+1 

i=l 

t (6-23) 

k 

6 dH r+1 = 6idH r + a°6 dHi + ^ j a^6 Q dH ni+ i, 

i=l 

where a s r = V r (n+k) ~ s .Of course formula (6.23) works separately for the geodesic Hamiltonians 
E r and potentials V r 



dE r+1 = N*dE r + c^dEx + afdE ni+1 , (6.24a) 

i=l 
k 

dV r+1 = L*dV r + a° r dVi + ^2 a^dV ni+1 . (6.24b) 

i=l 

Proof. The proof is inductive. Assuming that (6.23) and (6.24b) are valid for the k-hole case, 
with holes at positions ^( n + k )- n ^i = 1, we prove the validity of (6.23) for the (k + 1)- 
hole case with an extra hole at the position ^( n + k + 1 )- n k+i _ We ^ a ]^ e ^ e fc_hole and (k + l)-hole 
separability conditions as in (6.17) and (6.18). Our assumptions are as follows 

k 

dH r+1 = N*dH r + a^dHx + ^ a^dH ni+1 , (6.25a) 

i=i 

k 

da? +1 = L*dap + a°da? + ^ a T 3 ''dc£ +1 , (6.25b) 

3=1 

where i,j = 0, ...,k,no = 0. Notice that (6.25b) is a particular case of the condition (2.24) and 
follows from the fact that in a Stackel case otij are particular basic separable potentials. From 
(6.22) and (6.16) we find the following useful relations 

— ~ a ~ — 1 

H r+1 = H r q - — ff„ fc+1 _i, Hi = g F nfc+1 _i, r = 1, ...,n + k, (6.26a) 

a n k+1 -l a n k+1 -l 

K+f = op ~ r = 1, n + k, (6.26b) 

a n k+1 -l 

c%? = -P— , (6.26c) 

a n k+1 -l 
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a®, , i = a 



i V a" 1 a 



a 



t r+1 « r+1 -r ^ nk+i 

"n fc +i-l j =1 



r a? 



+1 ^0 " ^n k+ 

"fe+i-1 i=l 



-,-E 



(6.26d) 



where i7 ni = = 0, i = 1, a"; = 5^. Expressing by from (6.26a)-(6.26d) and 
using assumptions (6.25a), (6.25b) we find, after long but straightforward calculations, that 

k 

N*dH r + a° r dH l + Y,^ +1 dH nt+ 2 + a? k+1 dH nk+1+1 



i=i 



dH; 



H, 



ri fc+ i-l \"n fc+ i-l 



= d\H r - 



a n k+1 -l 



-H, 



n k+1 -l 



dH r +\ . 



Using the known relation for the Benenti potentials 

K r+1 = Y,p k L r -\ 

k=0 

one finds from (6.25b) that 

dV r+1 = [(L*Y + a\{L*) r - l + ... + a Q 1 ]dV l + [a^(L*) r - l + ... + a^} 
dV ni+1 + ... + K*(L*) r " 1 + ... + a^]dV nk+l 
= A®dVi + A^dV ni+1 + ... + A^dV nk+1 , 



where 



A~, — 



Pni Pni-1 " 

Pn k Pn k -1 

K r -1 



1 



Pm-i 



Pnu-l 



K r -k 
Pni—k 

Pn k —k 

■ Kr-k 
Prij—k 

Pn k —k 



(6.27) 



with the row (/£>„. _i, ...,p n ._f,) missing in A™* and again Ki in the calculation of determinants 
are treated as symbols not matrices. 

As in the one-hole case (5.21), the d p part of of (5.19a) gives us the analog of formulae (4.5) 
for the fc-hole case 



K r+1 = LK r + ap + 



i=i 



Then, from (6.28) we find that for L function (5.22) the following relation holds 

k 

L = E 2 - ot{Ei - ^ a^E ni+1 , 
i=i 



(6.28) 



(6.29) 
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hence 



k 
i=l 

k 

= + (6.30) 

i=l 

so obviously, L is not a conformal Killing tensor for G given by (6.8). 

Notice, that although the number k can be arbitrary large, nevertheless, the maximal number 
of nonvanishing terms d^ i dH ni J r \ in (6.23) is lower or equal to n. In fact, if rii and rij+i are two 
successive numbers, i.e. n« + i = + 1, then a^dH^+i = a^ l dH rH+1 = 0, as from construction 
H ni+1 = 0. Hence, for a string of successive numbers rii — Sj + 1, rii — Sj + 2, rij, only the term 
with nj is nonzero in formula (6.23), as 

H ni - Sl +i = H ni - Si+ 2 = ... = H ni = 0. 

Thus, assume that in the sequence 

jj^(n+k)-l + ^^(-+fe)-2 + + x n+k 

we have I strings of holes, where the i-th string has Sj holes and si + ... + si = k. Then, the 
quasi-bi-Hamiltonian chain (6.23) takes the form 

dH r+1 = N*dH r + a"Wi + <W ni+ i + ... + aP'dflnj+i, (6.31) 

where no = 0. 



6.6 Gel'fand-Zakharevich bi-Hamiltonian representation 



Now we lift our quasi-bi-Hamiltonian representation, constructed in a previous subsection, into 
a GZ form. Adding l + l Casimir coordinates q (with respect to 9q) and extending Hamiltonians 
H r to the form affine in Cj 

l+i 

h r = H r + Y^ a^Ci, (6.32) 
i=i 

one can transform the quasi-bi-Hamiltonian chain on T*Q into a bi-Hamiltonian chain on 
T*Q x R x ... x R, being a composition of I + 1 bi-Hamiltonian sub-chains 

s v y 

l+l 

TTodh ni = 
w dh ni+ i = TTidh ni 

: (6-33) 

Kodh ni+1 - Si+1 = 7Tl(^n i+ i-s i+ i-l 

= 7rid/i ni+1 _ Si+1 , 
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where i = 0, 1, I, h = ci, h nj = Cj+i, /i n;+ i-s !+1 = ^n ;+1 = hn+fc and 






... 
























/ 





/ 


01 


9odhi dh ni+ i . 


. 6odh ni +i 






-(9 dh 1 ) T 






7Tl = 




— (9odh ni+ i) T 











— (6odh ni+ i) T 







J 

Again 7To and 5fi are compatible Poisson structures as, according to Lemma 7, tt\ takes the form 

l+i 

TTl = TTlD + ^ Xj A Zi, 
i=i 

where Zi = -§^,Xi = TTodh ni+ i, Lz^q = 0. The reduction of / + 1 chains onto the symplectic 
leaf ci = ... = q + i = of 7To along the distribution Z = (^-, d ^ ) reconstructs immediately 
the quasi-bi-Hamiltonian chain (6.23). 

If we introduce the Poisson pencil tt^ = tt\ — Sttq, I + 1 chains (6.33) can be written in a 
compact form 

7^(0 = 0, h(0 = (hW(o,h (2) (0,---,^ l+1) (0) T , 



where 

h^(0 = h ni C i+1 - ni - Si+1 + ... + hn i+1 - Si+1 , i = 1,2, ...,/ + 1. 
The quasi-bi-Hamiltonian chain (5.17) takes the form 

e(dH(S)+a(0 o H w {q,p)=0 
t 

(TV* - HI)dH(0 + a(0dH (1) (q,p) = 0, 



(6.34) 
(6.34a) 



(6.35) 



where 



— Tlj — Sj+1 + 1 



+ ... + H n 



■s i+ i 



i = 1,2,. ..,1 + 1, 



H{i) — (Hi,H m+1 , ...,H ni+1 ) q 



and 



a(0j= g c . »,J = l,...,i + l. 



(6.35a) 



(6.35b) 
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The Stiickel separation conditions for extended Hamiltonians h r take the form (3.18b) 

Cl £(»+*) + h l ^ n+k ^ 1 + ... + V_ sl _ 1 ^ n+fc )-( ni - S1 - 1 ) 
+c 2 & +k) ~ ni + V+ie (n+fcK(ni+1) + ... + (6.36) 



i=l 

Remark 19 Systems considered in this paper, although obtained through the deformation proce- 
dure on the level of Hamiltonian functions, are far from being trivial generalizations of Benenti 
systems. There is no obvious relations between solutions of a given Benenti system and all its 
deformations. In each case we have a different inverse Jacobi problem to solve. Notice, that 
the common feature of appropriate deformed systems is the same set of separated coordinates, 
determined by the related Benenti system. 

7 Examples 

7.1 Henon-Heiles system 

Let us consider the integrable case of the Henon-Heiles system 

(q 2 h = ~qV, (7.1) 
with the corresponding Lagrangian 



L = \{q l ? t + \(q 2 ? t ~ {q l f ~ + cq\ (7.1a) 



The bi-Hamiltonian chain is of the following form [42] 

nodho = 

irodhi = X\ = iridha 
nodh2 = X2 = nidhi 
= v:\dh2 



(7.2) 



where 



h = c, 
1 

= 2 1 

= E 1 + V 1 (q) + Pl (q) c = H l +p l (q)c, 

h2 = ^ 2 -^ + ^( 9 2 ) 4 + |tf) 2 tf) 2 -^ 2 ) 2 
E 2 + V 2 (q) + p 2 (q) c = H 2 + p 2 (q)c, 



h = Ipl + lpt + iq'f + lq'iqr-cq 1 
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7T0 



#0 












#0 



/ 








1 


















1 






-i 













V 





-1 








/ 



7Tl 



01 



"l 


rlh \ 
UQUIll \ 


— (ovarii) 


n / 










W 


-I 1 ~h 


2 


-W o 


~\V2 



\a 2 \ 



\P2 



Notice that 

The quasi-bi-Hamiltonian chain is 



G 



L 



I 1 w 



dh r 



\a 2 

W -q 1 



dH r+1 = N*dH r + p-dHx, p r = — L , r = 1, 2, 

oc 

where N* = O^Oi. 

The transformation to separated coordinates (A,^) takes the form 



Pi 



A 1 - A 2 ' A 2 -A 



+ 



q 1 = X 1 + A 2 , g 2 = 2v C A T A 2 , 
A 2 /z 2 



l - 



P2 



Pi + P2 



A 1_ A 2 ' A 2_ A 1 



1 / ' 



and the separability conditions for Hi respectively hi are reconstructed from separation 

H l i+H 2 = l -^ 2 +e, 
c e + tut +h 2 = hp 2 + e. 



l-hole deformation of the Henon-Heiles. 

The only possibility is ri\ = 2. Then, 



G = G 

Pi 



4 o 



,K 2 = 0,K 3 = -Ki 



W?-{q l ? 



and 



Vi = --Vi, V 2 = 0, V 3 = V 2 -^V 1 , 
Pi Pi 



«r = Pr~ Pr-lP2Pl . Pr = Pr-lPl ■ 

The quasi-bi-Hamiltonian chain takes the form 

dH r+1 = N*dH r . + a r dH x + P r dH ni+1 , r = 1,2, 3. 
while two bi-Hamiltonian sub-chains are 



TTodho = 

TTodhi = X\ = TTidho 
= ix\dh\ 



7Todh 2 = 

TT dh 3 = X 3 = TTidh 2 
= TT\dh 3 
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where 



ho 
hi 

h 2 
h 3 



ci, 

5^ + 5^ + ta') , + 5tf)*-c[B + ^tfffl-4 

Ei+VqlW+a^q) c 1 +P 1 (q)o 2 , 
c 2 , 



1 2 1 <?2 2 1 

E 3 + V 3 (q) + a 3 (q) a + (3 3 (q)c 2 



1 (g 2 ) 4 1 (g 2 ^ 2 
16 4 q 1 



and 



#0 


\ 


















Oodhi 9odh 3 \ 







Til) - ■ 7T| - ifl„f///| )' 

V -(MM T 

The transformation to separated coordinates is given by (7.4) and the separability conditions 
for Hi respectively hi are represented by separation curves 



H l e+H 3 = -^ 2 +i\ 



C1 e + h ± e + C2 i +h 3 = + 1. 



The related Euler-Lagrange equations (3.2) arc 

iWu + ~ \i<£? = + [1 - \{ q l )-\ q 2 ?]c. ~ ( 9 1 )- 2 c 2 , 



1 



Q ' Qu + tilt = ~Q + g(9 )" Vci- 

2-hole deformation of the Henon-Heiles. 

The only possibility is 2-hole string m = 2, n 2 = n\ + 1 = 3. Then, 



(7.8a) 
(7.8b) 

(7.9) 



G 



P1-P2 



-G 



W7W 





4 



K 2 = K 3 = 0, 



Ka = Kn 



■\q\q 2 f i(g2)S + i (9 i)2 + g2 



1\3 



Fi = -^— Fi, y 2 = y 3 = o, y 4 = y 2 -^ 2 -y 1 , 



Pi - P2 



Pi - P2 



p 2 Pl 

a.\ = p 1 5 , a 2 = a 3 = 0, a 4 — 9 

Pi - P2 



p\ 



01 — 2 j 

Pi- to 



Pi ~ P2 

p 2 = p 3 = o, p 4 = J^±- 

Pl - P2 
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and two bi-Hamiltonian chain are 

TTQdho = 



nodhi = X\ = iTidho 
= TT\dhi 



TTodhs = 

TTQdh4 = X4 = TTldhs 

= TTidhi 



(7.10) 



where 



h 
hi 

h 2 
h± 



ci, 



4(g 1 ) 2 + q 2 



[2p( + 2p I 2 + A{q'Y + 2q'( q I Y}-ciqi 



4c 2 



4(g 1 ) 2 +q 2 ' 



0, 

C2, 

4(gl) 2 + g2 [^V)M + PfoW + ^ 2 ) 3 )P1P2 " (?V) 2 



+2(g 1 ) 3 )^ + 1(^)6 



1 

16 



ci(q 2 ) 6 + c 2 q\q'f}, 



l ( '„2\2i 



7T0 



#0 


\ 









J 










( 


Oi 


9 dhi 9 dfi4 


TTl = 




-%dhi) T 







\ 


-(Oodh 4 ) T 



The transformation to separated coordinates is given by (7.4) and the separation curves for Hi 
respectively hi are 

(7.11a) 



cit + h^ + ^i + h^^ + e- 



(7.11b) 



7.2 7-th order stationary KdV 

Let us consider the so-called first Newton representation of the seventh-order stationary flow of 
the KdV hierarchy [43], [22]. It is a Lagrangian system of second order Newton equations 

4 = -io(^) 2 + V 

it = -l6q 1 q 2 + W(q 1 f + 4q 3 (7.12) 
q\ = -20gY - 8(q 2 ) 2 + 30(g 1 ) 2 g 2 - ^(q 1 ) 4 + c 
with the corresponding Lagrangian 

£ = qkt + Wtf + 4 <?V " lO^V - Sq\q 2 ) 2 + lU(q l fq 2 - S^ 1 ) 5 + cq\ (7.12a) 

so that it can be cast in a Hamiltonian form. In fact, the above system can be represented as a 
bi-Hamiltonian vector field Xi belonging to a bi-Hamiltonian chain of the form 

nodho = 

wodhi = Xi = nidho 

wodh2 = X 2 = nidhi (7.13) 
nodh 3 = X 3 = TTidh 2 
= iridh? 
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hi 



P1P3 + \p\ + lO^V - Vg 3 + Sq\q 2 ) 2 



10(0 V + 3(9 1 ) 5 -at 
Ei + Vi(q) + Pl (q) c = Hi+ Pl (q)c 

kVs - Wvi + h^ps - biP2 - tfpm + V) 2 + kov 

§ (g 1 ) 6 - 2(g 2 ) 3 + (, 3 ) 2 - 6, W + Ictf) 2 + icg 2 
E 2 + V 2 {q) + p 2 {q)c = H 2 + p 2 {q)c 

1{Q 2 ) 2 P 2 3 + W?Pl + hi + \£v\Pi + \q 2 P\P3 - h\ 2 P2Pz 



h 3 P2P3 - 3(0 V) 2 + qHq 2 f + ftf) V + WW 
Kg 1 ) V + (q 2 ) V " (^W - Wq 2 - \a? 



= E 3 + K 3 (g) + P 3 (?)c = H 3 + p 3 (g)c, 
with the corresponding operators 7To and tti 

n n = 



n 



1 — o 



#0 


06x1 


•)• 


o = 


( 


3 


0lx6 









-h 











q 1 


-1 














q 2 





-1 











2q 3 


q 2 


q 1 


-q 1 


-q 2 


-2g 3 





P2 


P3 


1 





-q 2 


-P2 











1 


-q 1 


~P3 












-2X 


T 
i 







h 



2X, 



and Xi = Oodhi, d = 



d_ d_ 

dq ' dp 

Hamiltonian chain is of the form 



/ 



so that the operator 61 in the corresponding quasi-bi- 



/ 











q 1 


-1 





\ 













q 2 





-1 















2q 3 


q 2 


q 1 






-q 1 


-q 2 


-2q 3 





P2 


P3 






1 





~q 2 


~P2 










V 





1 


-q 1 


~P3 








/ 



From the form of hi one can directly see that the inverse metric tensor G expressed in (q,p) 
variables has in this example an anti-diagonal form 



G 
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while the conformal Killing tensor L has the form 



and hence K\ = I, A\ = G, 



*-3 



/ q 1 -1 \ 
q 2 -1 

-.3 „2 „1 



\ 2q 6 q 1 





1 

2 


( -q 1 


-1 





K 2 = 


q 2 


-2q x 


-1 








? 2 








/ o 


-1 




A 2 = 


1 

2 


-1 


-2^ 


q 2 






W 


q 2 


2q 3 



-q 1 q 2 - 2q 3 
(q 2 ) 2 




1 



1 



1 ! T 1 ?/ 



g 2 

g 1 *? 2 — 2q 3 



q 1 

{q l ? 

q 2 -q 1 q 2 - 2q 3 (q 2 ) 2 J 

The quasi-bi-Hamiltonian chain is given by (7.3) with r = 1,2,3, where 

pi = -q\ p2 = \(q 1 ) 2 + \q 2 , p s = -\q l q 2 -\q 3 - 

The transformation (A,//) — > (<?,p) is constructed from the relations 

(? l = A l + A 2 + A 3 ) l ((? l ) 2 + l g 2 = A l A 2 + A l A 3 + A 2 A 3 ) ^y + ^3 = ^2^ 

(the explicit formulas are given in [22]) and the separation curves for Hi respectively hi are 



Hit 2 + H 2 Z + H 3 = -p 2 + 16£ 7 , 
c£ 3 + h^ 2 + h 2 £ + h 3 = In 2 + 16£ 7 . 



(7.14a) 
(7.14b) 



1-hole deformation. 

There are two admissible cases of one-hole deformation: n\ = 2, 3. Here we illustrate the 
case n\ = 3. The deformed Hamiltonians are the following 

1 1 T ~ 

= -p T Gp+V 1 , 









P2 




= #1- 


% 


= 0, 


H A 


= H 3 - 



Pi 

P2 



1 



(7.15) 



where 



G = — —K 2 G 

P2 



^H 2 = V, 

P2 2 



1 



1 



1 u 



z 1 \ 

„2 



(g!) 2 + 2g 2 | ; 2 
9 -9 



Ai = G, 



1 ~ 2 -k 3 / 
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Ao 



--{K 2 - K x Kz)G 
92 

( 



(g 1 ) 2 + 2q 2 




2q 1 



2q 1 (q 1 ) 2 - 2q 2 \ 

3( (/ 1 ) 2 - 2q 2 -2q l q 2 



\ (q 1 ) 2 - 2q 2 -2q 1 q 2 -4q 1 q 3 / 



^3 

A, 



= 0, 



= —KZG=- 



1 1 
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4 (q 1 ) 2 + 2q 2 



( (q 1 ) 2 + 2q 2 {q 1 ) 3 - 2q 3 

* q 1 ^ 1 ) 3 -2q 1 q 2 - 4q 3 } 

* * 



2(g 2 ) 2 - {q l ?q 2 ~ 2q 1 q 3 
-(q^q 2 -2{q 1 ) 2 q 3 -2q 2 q 3 



\ 



2(q 2 ) 3 + (q 1 ) 2 (q 2 ) 2 +4q 1 q 2 q 3 +4(q 3 ) 2 ) 

the respective potentials are 

S^ 1 ) 6 - 8(q r ) 2 (q 2 ) 2 - 10(qyq 2 + 8(q 2 ) 3 - 4(q 3 ) 2 + 24q 1 q 2 q 3 

v 2 

% 



(q 1 ) 2 + 2q 2 

[(g 1 ) 7 - S^JV " BfoW + 2(g 1 ) 3 (g 2 ) 2 + 4(g 1 ) 2 gV ~ *Z V) 3 
-2q\q 3 ) 2 + 4{q 2 ) 2 q 3 ]/U) 2 + 2q% 
0, 

wftf? - ^ i ne? + 2( 9 i ) 3 (g 3 ) 2 - gVfo 3 ) 2 + ^ w 

-5(g 1 ) 2 (g 2 )V + (g 3 ) 3 ]/[(g 1 ) 2 + 2g 2 ], 

and (*) are elements making matrix symmetric. The quasi-bi-Hamiltonian chain takes the form 
(7.6) with r = 1, ...,4 and 

_ _ 9r-lPz n _ Pr-1 

oi r — p r , p r — , 

92 92 

while the related bi-Hamiltonian chain splits onto two sub-chains of the form 



TTodh3 = 

= TT\dhi , 



where 



TTodho = 

nodhi = X\ = TTidho 
TTodh2 = X2 = TTldhl 
= TT±dh2 

1 7 

^o = ci, /13 = c 2 , /ii = -p A, p + Vi + OLiCi +/3jC 2 , z = 1,2,4. 



(7.16) 



Separation curves for respective Hamiltonian functions are 

1 



#i£ 3 + H 2 e + H 4 = -fi 2 + 16£ 7 , 
cit + h^ 3 + h 2 f + c 2 t + h A = Ifi 2 + 16£ 7 . 
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2-hole deformation. 

There are three admissible cases of 2-hole deformation. One related to 3-Casimir Poisson 
pencil: n\ = 2, n 2 = 4 and two cases related to 2-Casimir Poisson pencils, i.e. n\ = 2, n 2 = 3 
and n\ = 3,n 2 = 4. Here we present the first case. The deformed Hamiltonians are the following 



Hi 

H 2 
H 3 

H A 
H 5 



H 2 = \ p T Gp + V u 

PlP2 ~ P 3 2 



= 0, 



Hi + P2 Pi ^ = 1 T Aa p + y^ 
PlP2 - P3 2 

0, 

#3 + 



PlP2 ~ P3 



where 



A x = G = 

A 2 = 0, A 3 
Vi = 



PlP2 - P3 
1 

PlP2 ~ P.! 

4 



-K X K 2 G = 



Qi ~ QiQ2 + 93 



/ 1 
1 



V q 1 -g 2 - 




K 2 (K 3 - Kl)G, A 4 = 0, A 5 



PlP2 ~ P 3 



-K 2 KiG, 



q3 - 9i92 - 9i 



2(7 2 - 3(/? 

3V2, V 2 = 0, V 3 = Vi + ^ ^Va, 



93 - 9i92 - 9i 



^4 = 0, V 5 = V 3 + ^±4 y 2 . 

93 - 9192 - 9l 

The quasi-bi-Hamiltonian chain takes the form 

ciff r+ i = iV*dF r + a°(iFi + a^d^ 3 + ajdH 5 , r = 1, 5, 



(7.17) 



(7.18) 



aj; = p r + 



Pr-2P2P3 ~ Pr-lP2 



Pr-lP2 ~ Pr-2P3 



P\P2 ~ P3 PlP2 ~ P3 

and the respective bi-Hamiltonian chain splits into three sub-chains 



Pr-2P\ — Pr-l 
PlP2 - P3 



nodho = vr d/i2 = 7r d/i4 = 

TT dhi = X\ = TTidliQ 7Todh 3 = X 3 = TTidh 2 ttq^s = X$ = TX\dh 4 



= ix\dh\ 



where 

























and 



h = ci , h 2 = c 2 h 4 



= TT\dh 3 



= -K\dh^^ 



7Tl 



01 


Oodhi 9odh 3 Oodh^ 


-(MM T 




-{0 dh 3 ) T 





-(6 dh 5 ) T 





C3, /tj = #j + q°ci + afc2 + ajc 3 , i = l,3, 5. 



(7.19) 
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The respective separation curves are 

Hi^ + + H 5 = ^ 2 + 16£ 7 , (7.20a) 

Cl £ 5 + + c 2 f + h 3 e + c 3 £ + h 5 = l -tx 2 + 16£ 7 . (7.20b) 

8 Summary 

We have presented a geometric separability theory of Liouville integrable systems with n quadratic 
in momenta constants of motion 

Hi(q,p) =p T Mq)p + Vi(q), i = l,...,n (8.1) 

and with separation curves of polynomial type 

HiC 1 + - + HnC" = \f (0^+^(0, m n = < m n _i < ... < mi £ N, (8.2) 

where /(£), and 7(£) are Laurent polynomials of £. Such systems can be lifted to the class 
of Gel'fand-Zakharevich bi-Hamiltonian systems, defined by linear Poisson pencils and their 
Casimirs: polynomial functions in pencil parameter. First we briefly summarized the geometric 
separability theory of all GZ bi-Hamiltonian systems. Then, we reviewed with details a special 
class of such systems introduced by Benenti and characterized by separation curves of the form 
(8.2) with rrik = n — k. Finally, we proved in a very systematic way that all remaining separable 
systems, characterized by separation curves of the form (8.2), are constructed by the appropriate 
deformations of related Benenti systems. The most fundamental and surprising result of the 
presented paper can be formulated in the following way. If, for a given canonical coordinate 
system (qi,Pi),i = l,...,n, we have a pair of objects, i.e. contravariant metric tensor 
and related special conformal Killing tensor L, then we can construct systematically, in these 
coordinates, all Liouville integrable and separable Hamiltonian systems (8.1) with respective 
separation curves of the form (8.2) and explicit form of transformation to separated coordinates. 
Observe, that accordind to what was said in Section 4, the separation curve of geodesic motion 
for Gr(°) is the following 

EiC' 1 + E 2 C- 2 + - + E n = i/i 2 . (8.3) 

So, the passage from the system (8.3) to the systems (8.2) is constructive and determined 
completely by the pair (G^°\ L). 
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